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Abstract 



^ ■ Presented is the exact solution of the problem of finding a potential of an infiaton scalar 

^ . field for which adiabatic perturbations generated during a de Sitter (infiationary) stage in 

the early Universe have the exactly fiat (or, the Harrison- Zeldovich) initial spectrum. This 
solution lies outside the scope of the slow-roll approximation and higher-order corrections to 
00 ' it. The potential found depends on two arbitrary physical constants, one of those determines 

an amplitude of perturbations. For small (zero) values of the other constant, a long (infinite) 
. infiationary stage with slow rolling of the infiaton field exists. 

; PACS numbers: 98.80. Cq, 98.80.Hw 

^ ■ Great increase of amount and accuracy of cosmological observational data obtained al- 



O 

o 



iT) . ready and expected in near future makes real the ambitious programme of determination 

of an initial power spectrum Po{k) of inhomogeneous density perturbations in the Universe 
directly from data. Then, if working within the scope of the simplest versions of the infiation- 
ary scenario of the early Universe (i.e., with one slowly rolling scalar field (p - the infiaton), a 
natural next step is to reconstruct an effective self-interaction inflaton potential V{(f)) lead- 
^ . ing to the generation of such Po{k) from inflaton quantum vacuum fluctuations during a de 

Sitter (inflationary) stage. Note that in inflationary cosmology, an "initial" power spectrum 
means the spectrum that had been formed by the end of the inflationary stage, i.e. by the 
^ ' beginning of a post-inflationary power-law evolution of the Universe (it was flrst calculated 

■ in [1] for the R + E? model [2] and in [3, 4, 5] for the "new" inflationary model [6, 7]). 

Here k = |k|, and the spatial dependence exp(ikr) of all Fourier modes of adiabatic (scalar) 
perturbations is assumed. 

From the mathematical point of view, since quantum generation of perturbations re- 
duces to a kind of scattering problem, the reconstruction of V{(f)) may be considered as 
a speciflc inverse scattering problem. Its solution is known in the slow-roll approximation 
for the inflaton fleld, including flrst [8], second [9] and third [10] order corrections to this 
basic approximation with respect to small slow-roll parameters, and in the so-called general 
slow-roll approximation [11, 12]. It is non-unique if additional information about the spec- 
trum of primordial gravitational waves generated during a de Sitter (inflationary) stage [13] 
is not used. Namely, generally there exists a one-parameter family of V(0) producing the 
same Po(^) (here and below I don't count one more trivial free parameter corresponding to 
invariance with respect to a constant shift of the inflaton fleld: 0^0 + 0o)- However, an 
exact solution of the reconstruction problem is desirable for at least two purposes: 

1) to investigate how good is the convergence of perturbation series in powers of slow-roll 
parameters (if it takes place at all); 

2) to determine the exact degree of degeneracy of the problem, i.e. to flnd the measure of a 
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set of potentials producing the same perturbation spectrum. 

In particular, the problem of accuracy of the slow-roll approximation prediction for Po{k) 
(including higher order corrections) has been intensively and critically studied recently us- 
ing different methods: [14], [15] (the uniform approximation), [16] (the improved WKB- 
approximation) and others. 

By an exact solution I mean a solution of the following system of equations for a spatially- 
flat Friedmann-Robertson- Walker (FRW) background with a scale factor a{t) and scalar 
(adiabatic) perturbations described by the Mukhanov variable Q = u/a: 

(1) 




^ + 3ff^,+ ^ = 0. (2) 

drf \ z drf ) ' 
obtained without any approximations. Here 

d a0 r dt 

H ^ - , — , r]^ —- , 4 

a H J a[t) 

dot means the derivative with respect to t, Uk{ri) exp(ikr) is the wave function of a Fourier 
mode of the quantum field u (the c-number multiplying the Fock annihilation operator dk), 
and c = h = \ IS put throughout the paper. The variable Q [17] is equal to 5(j)L + in the 

longitudinal gauge ($ is the quasi-Newtonian gravitational potential), or to 505 — g^(/x-|-A) in 
the synchronous gauge (/x and A are the Lifshits variables). The normalized initial condition 
for Uk corresponding to the adiabatic vacuum at i — > — oo (77 — > —00) is 



At late times during an inflationary stage in the super-horizon regime {k <^ aH, r] ^ 0), 



^ const = ({k) (6) 

(C = —h/2 in the notation of [4]). 

Then the initial spectrum of adiabatic perturbations for a post-inflationary cosmology in 
the super- horizon regime is (assuming the absence of non-diagonal pressure components): 

Here t = corresponds to the end of inflation. For historical reasons, the slope of the 
spectrum is deflned with respect to density perturbations in the non-relativistic dark matter 
-I- baryon component at the present time, {5p)k ~ —k'^^k/^^nGa'^ before integration over d^k. 
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So, ns = 1 + "^'^y . Finally, using the equation H = -AnGcj? that follows form Eqs. (2) 
and (3), Eq. (2) can be recast in the Hamilton- Jacobi form [18]: 

where the prime denotes the derivative with respect to 0. 

Exact solutions of the inverse problem of reconstruction of given Po(^) are known for 
the following two cases only, if not speaking about solutions describing universes collapsing 
towards a singularity. 

1) A power-law perturbation spectrum with the slope ng — const < 1 [19]. Then 



V{(f)) oc oc exp ( ± ) , a{t) oc ^^ g = > 1 . (9) 

V V / i- - ns 

This is just the power-law inflation. Considered as a function of (pit), H is related to 
V{4>) through Eq. (8). Note, however, that this is not the only potential producing the 
ns — const < 1 spectrum. 

2) The case when no perturbations are generated at all (no real created quanta of the inflaton 
field) [20]: 

i/(0) = H, exp(27rG02) , V(0) = ^ (l - exp(47rG02) . (10) 

In literature, this case is sometimes incorrectly referred as the potential generating the 
^5 = 3 perturbation spectrum. However, one should not forget that generated perturbations 
are quantum (even quantum-gravitational) and require renormalization. After subtraction 
of the vacuum energy uj{t)/2 = k/2a{t) of each mode, no created fluctuations remain in 
this case. Moreover, a number of real inflaton quanta generated in each perturbation mode 
k should be large, because in the opposite case they may not be interpreted as classical 
perturbations after the end of inflation (see [21] for a more detailed discussion of this point). 

Strictly speaking, there is no exit from inflation for the potential (9), and the potential 
(10) does not admit a low curvature regime at all. However, in the former case V{(j)) can be 
deformed such that it reaches zero at a sufficiently large value of 4>. This will result in a very 
small change of the perturbation spectrum at present scales of interest that may be safely 
neglected. Sometimes, the case of a parabolic potential near its maximum ^(0) — Vq — ^^^^ 
is mentioned as an exactly soluble case. However, it is not such the one in our terminology 



since in this case H{(f)) is approximated by the constant value Hq = ySnGVo/S. 

In this paper, a family of exact solutions for the case = 1 is constructed. It is just the 
initial spectrum proposed by Harrison and Zeldovich [22] , after all, for beauty reasons. Note 
that it satisfies the most recent CMB data [23, 24]. Let us first consider what follows for 
this case from the slow-roll approximation. Then, the leading term in the power spectrum 
reads 

k'Cik) (X (^^^ ^ ^ , (11) 

where tk is the moment when k = aH. It is clear that, to get ns = 1, V^^^/V should 
not depend on 0. Therefore, V{4>) oc Note that this solution of the reconstruction 
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problem is unique for a given amplitude of the flat spectrum. This kind of inflation was 
dubbed intermediate inflation in [25] (see also [26]). Its scale factor behaviour is a{t) oc 
exp (const ■ nee more, it does not have an exit from inflation, so it should be modified 

at large (I). A next order slow-roll correction to this potential was considered in [27]. 
To obtain an exact solution for i/(0) and V{(f)) in the case ns = 1, note first that, for 

-t4 = 4' (12) 

Eq. (3) reduces to the equation for a massless scalar field in the de Sitter background and 
has the solution 



satisfying the initial condition (5). Let us write the general solution of Eq. (12) in the form 

z = -^fl + ^V r/<0, (14) 

where A, rjo are constants. The limiting case rjo 0, when the first term in brackets may 
be neglected, is not interesting because it corresponds to a collapsing universe (however, 
it is "dual" to the case 770-^00 considered below). The power spectrum of the growing 
perturbation mode is Po{k) = l/An'^B'^ and does not depend on rjo {rjo appears in the 
amplitude of the decaying mode only and makes it non-scale-free). Thus, we have got the 
exactly flat spectrum. Present observational CMB data [24] fix the quantity B with 10% 
accuracy: 

where A is the quantity introduced in [24] and r is the optical length after recombination. 
In this notation, A = 0.9 corresponds to the value A = 4.3 • 10~^ of the other quantity A 
introduced in [28] to characterize an amplitude of initial perturbations (and conjectured to 
he in the range (3 — 10) • 10~^ in that paper). 

Since the aim of this paper is to find some exact solution, I will not investigate if there 
exist other forms of z leading to the ns = 1 spectrum, too. The absence of other solutions 
for z would immediately follow from scaling arguments if we assume that Uk oc k^^^'^f{kr]) 
for all 77. However, the latter assumption might not be necessary. Moreover, I will consider 
only one particular case of Eq. (14) corresponding to the limit 770 — > 00. 

So, let z — —B/t). Let us express all quantities of interest as functions of 0: 



= -4nG ^ , Ina = I H{t) dt = -4nG j ^ d(f) , 



(16) 
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Equating the last line in Eq. (16) to —B/r), we get the following equation: 

I P{<t>) dcP = -BHP , P = 1^ exp (attG J§jd(l>^ . (17) 
After differentiation, Eq. (17) reduces to P = -B{HP' + H'P), or 

Let us introduce dimensionless variables 

X = V4^0 , y = bV4^H , v{x) = V{(j)) . (19) 

Then, from (8), v ^ - {l/3){dy/dx)^. For these variables, Eq. (18) reads: 

d^y / dy\^ ^ dy ^ 2 /2q\ 
dx'^ \dx J dx 

After dividing by y"^, the last equation can be integrated to dy/dx — xy — 1 (an integration 
constant is excluded by shifting x, i.e., 0). Therefore, 

y ^ e^y^ ^l^-^y^x + , (21) 

where C is another integration constant. This just yields us a one-parameter family of 
solutions having ^5 = 1. The so-called slow- roll parameters for this solution: 

AttG H y dx^ y 

The partial solution with C — has an infinite inflationary stage which is just described 
by the slow-roll approximation for x ^ 1. Its graph is plotted in Fig.l. Its large-x expansion 

1 1 3 15 17 9 

•X/ JU JU JU KJJU JU 

It is straightforward to check that it leads to ^5 = 1 (as it should be) for the first [8] and 
second [9] order corrections to the slow-roll approximation. However, these corrections miss 
the whole 1-parametric family with C 7^ completely. 

For a; < 0, the solution with C = has rather peculiar behaviour: the potential v{x) 
reaches the maximum value Vmax ~ 7.252 at x ~ —1.326, becomes zero at a; ~ —1.618 and 
then going to —00 at x — > —00 (however, such effective potentials are considered in string 
inspired models now). In the latter limit, y — > 00, so we get an initial curvature singularity 
at a finite proper time to < 0. If t = is the moment when a; = {v{fd) = f — | ~ 1.237) and 
the infiationary stage begins, then |to| ~ i?~^(0) ~ BG^^^. The scale factor reaches zero very 
slowly: a{t) oc | ln(t — to)!""*^^^ for t — > Iq. Still the Riemann tensor is not twice integrable for 
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Figure 1: The dimensionless potential v{x) for the C = case. 



t —>■ to, so this singularity is a strong one. The same refers to all initially expanding (y > 0) 
solutions with C 7^ and C > — -\/27r - they all begin from such a singularity. 

By taking C < and very small, it becomes possible to construct a solution with a long 
but finite inflationary stage. Namely, if C = —\^x~['^exp{—Xi/2) with xi ^ 1, then v{x) 
becomes zero at x = Xi {y still remains ~ Xi^). In this case inflation ends (e, \fi\ ~ 1) at 
X = xi — 0{xi^). The total number of e-folds is A^tot = 27iG(f)1 = x\/2. Thus, C ~ exp(— iVtot) 
that is in agreement with the general principle that terms not caught by an arbitrary order 
of a WKB-type expansion are exponentially small. For x > Xi, one may put v = Then 
the kinetic dominated phase a{t) oc t^^^ follows the inflationary stage. Or, we may assume 
that V has a local minimum w = | — Xi)'^ around this point. It results in oscillations in 
(j) and the matter-dominated post-inflationary stage a{t) oc t"^/^ . 

Finally, note that the spectrum of gravitational waves (GW) is not flat for this model: 
for 1 <^ X <^ xi, the tensor-scalar ratio and the slope of the GW initial power spectrum 
r = — 8nT = = d>/N where N is the number of e-folds from the beginning of inflation. 

The present upper observational bound r < 0.36 [29] requires N > 22 for the comoving scale 
crossing the Hubble radius at present. So, Ntot should exceed ~ 70 in this model. 

The research was partially supported by the Russian Foundation for Fundamental Re- 
search, grant No. 05-02-17450, by the Research Programme "Elementary Particles" of the 
Russian Academy of Sciences and by the scientiflc school grant No. 2338.2003.2 of the 
Russian Ministry of Education and Science. 
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